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Abstract. Recently it was shown that the excess of diffuse Galactic gamma rays above 1 GeV 
could be interpreted as a Dark Matter annihilation signal. From the spectral shape of the excess 
it is possible to determine a range for the allowed WIMP mass which can be used to test the 
supersymmetric parameter space. It is found that the EGRET excess combined with electroweak 
and other constraints, e.g. relic denisity or direct detection limits, is fully consistent with the 
minimal mSUGRA model for scalars in the TeV range and gauginos below 500 GeV. 

PACS. 12.60.Jv Supersymmetric models - 95.35.+d Dark matter 



1 Introduction 

Cold Dark Matter (CDM) makes up 23% of the en- 
ergy density of the universe [1] . One of the most pop- 
ular CDM candidates is the neutralino, a stable neu- 
tral particle predicted by Supersymmetry. The neu- 
tralinos are spin 1/2 Majorana particles, which can 
annihilate into pairs of Standard Model (SM) parti- 
cles. A large fraction of the annihilations is expected 
to go into quark- antiquark pairs. Since the DM parti- 
cles are strongly non-relativistic, the produced parti- 
cles are mono-energetic. In a recent paper we showed 
that the observed excess of diffuse Galactic gamma 
rays above 1 GeV, measured by the EGRET experi- 
ment, has all the properties of the ir° decays of such 
mono-energetic quarks originating from the annihila- 
tion of neutralinos with a mass around 60 GeV [2] ■ It 
was also shown that a newly determined extragalactic 
gamma ray background including Dark Matter anni- 
hilation (DMA) shows a less significant bump at the 
same energies [3]. 

It is the purpose of the present paper to see if this 
intriguing hint of DMA is compatible with Supersym- 
metry. Here we will concentrate on the Minimal Su- 
persymmetric Model with supergravity inspired sym- 
metry breaking (mSUGRA model). We assume that 
the EGRET excess originates from the annihilation 
of the stable, neutral lightest supersymmetric parti- 
cles, the neutralinos. Their mass is then constrained to 
be between 50 and 100 GeV from the EGRET data, 
which strongly constrains the masses from all other 
SUSY particles, if mass unification at the GUT scale 
is assumed. The upper limit on the WIMP mass de- 
pends strongly on the background model of the Galac- 
tic gamma rays. If one uses a conventional galactic 
model which reproduces the locally measured electron 
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and proton fluxes as well as secondary to primary ra- 
tios, one finds an upper limit of ~ 70 GeV. If one uses 
a model which was optimized to fit the gamma ray 
excess, it is found that the excess is still existent but 
not as strong as before, which leads to a higher mass 
limit of ~ 100 GeV. It will be shown that combin- 
ing the EGRET data with other constraints, like the 
electroweak precision data, Higgs mass limits, chargino 
limits and relic density leads to a very constrained 
SUSY mass spectrum with light gauginos and heavy 
squarks and sleptons @]. 

2 Predictions from mSUGRA 

The mSUGRA model, i.e. the Minimal Supersymmet- 
ric Standard Model (MSSM) with supergravity inspired 
breaking terms, is characterised by only 5 parameters: 
m o, mi/2, tan/3, sign(/x), Aq. Here mo and m\/2 are 
the common masses for the gauginos and scalars at the 
GUT scale, which is determined by the unification of 
the gauge couplings. Gauge unification is still possible 
with the precisely measured couplings at LEP . The 
ratio of the vacuum expectation values of the neutral 
components of the two Higgs doublets in Supersym- 
metry is called tan (3 and Aq is the trilinear coupling 
at the GUT scale. We only consider the dominant tri- 
linear couplings of the third generation of quarks and 
leptons and assume also Aq to be unified at the GUT 
scale. Electroweak symmetry breaking fixes the scale 
of /it, so only its sign is a free parameter. We use the 
positive sign, as suggested by the small deviation of 
the muon anomalous moment. 

The dominant annihilation diagrams of the lightest 
supersymmetric particle (LSP) neutralino are shown 
in figure [TJ The cross sections are proportional to the 
final state fermion mass, which originates either from 
the Yukawa couplings for the Higgs exchange diagram 
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Fig. 1. The dominant annihilation diagrams for the light- 
est neutralino, which is a linear combination of the gaugino 
and Higgsino states: |xo) = iVi|Bo) + N 2 \W$) + NslHx) + 
Ni\H2)- The dependence of the amplitudes on masses and 
neutralino mixing parameter Ni has been indicated. 



or from the helicity suppression at the low energies 
involved in cold DMA. Therefore heavy fermion final 
states, i.e. third generation quarks and lcptons, are ex- 
pected to be dominant. The Higgs exchange diagram is 
in addition proportional to tan (3 for down type quarks 
and l/tan/3 for up type quarks, indicating that top 
quark final states are suppressed for large tan (3. 

The W- and Z-final states have usually a much 
smaller cross section due to the weak couplings in- 
volved and are in addition kincmatically suppressed for 
the 60 GeV neutralino mass preferred by the EGRET 
data. 

For the pseudoscalar Higgs and sfermion exchange 
diagrams the annihilation rate, which is proportional 
to the cross section multiplied by the relative neu- 
tralino velocity, is practically independent of the cen- 
ter of mass energy, but the other diagrams show a 
strong energy dependence @]. This implies that for 
the present temperature of the universe close to ab- 
solute zero the neutralino annihilation is dominated 
by either sfermion exchange or pseudoscalar Higgs ex- 
change. The sfermion exchange is suppressed for the 
following reason. The Born mass of the lightest Higgs 
is below the Z° mass, but radiative corrections can 
boost it up to 130 GeV in the minimal mSUGRA 
model. These corrections depend on the heavy par- 
ticles coupling to the lightest Higgs, like the top and 
stop quarks. These scalars have to be sufficiently heavy 
in order to reach a Higgs mass above 1 14 GeV, which is 
the present lower limit from the direct searches at LEP 
[5J. Note that this is the limit on the Standard Model 
Higgs particle, but for heavy scalars the lightest SUSY 
Higgs particle has very much the properties of the SM 
Higgs, so the above limit is also valid in our case. For 
light ncutralinos, i.e. small m^, the Higgs mass limit 
excludes scalar masses below the TeV range, as indi- 
cated in Fig. [51 by the almost vertical line, labeled rrih- 
In addition, the excluded regions from Br(B — > Xsj) 
and the anomalous magnetic moment of the muon have 
been indicated (left from the corresponding lines). The 
mass spectrum has been calculated with the Suspect 
program 10J while the quantities fih 2 , Br(B — > X s j) 
and Z\a?, USY were calculated with micrOMEGAs 7 : . For 
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Fig. 2. The shaded area (blue) indicates the 95% CL 
parameter range in the mo-m^-plane allowed by the 
EGRET data, if the constraints from electroweak data, a 
neutral LSP and electroweak symmetry breaking (EWSB) 
are imposed as well. The individual constraints have been 
indicated by the lines and dots. The values of Ao = and 
tan/3 = 52.2 were chosen to be compatible with WMAP 
for the parameters marked by the star. 



the exclusion limits the following inputs were used: a) 
Br(B -> X s j) = (3.43 ±0.36) ■ 10~ 4 , which is the aver- 
age from BaBar, CLEO and BELLE QI]. b) the devi- 
ation of the anomalous magnetic moment of the muon 
a M from the expected value in the Standard Model was 
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The constraints on the supersymmetric parameters 
space are practically independent of Aq due to a co- 
incidence from the constrains from the Br(B — > X s j) 
rate and the lower limit on the Higgs mass of 114 GeV 
[5] . The region of large mo , for which no electroweak 
symmetry breaking (EWSB) is possible, has been in- 
dicated in Fig. as well as the region of small mo, 
where the stau would be the lightest SUSY partice, 
which is excluded, since the DM candidate has to be 
neutral. The mo values in the TeV range between the 
Higgs mass limit and EWSB limit are allowed by all 
constraints considered sofar. 

An independent check that the sfermion exchange 
diagram is strongly suppressed comes directly from the 
EGRET data: if the scalars are light, the stau is usu- 
ally the lightest scalar, in which case the stau exchange 
in the t-channel (left diagram of Fig. [1]) would be dom- 
inant, thus leading to tau final states. The low de- 
cay multiplicity of tau leptons leads to a much harder 
gamma ray spectrum from the hadronic decays, which 
is excluded by the EGRET data. 

For the allowed region in Fig. [2] discussed sofar the 
relic density constraint has not yet been considered 
and the interesting question is: does this rather nar- 
row region yield the correct relic density? The answer 
is simple: almost any pair of mo and mi/2 values are 
allowed by the relic density, if the values of the re- 
maining parameters, namely Aq and tan/3 are chosen 
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Fig. 3. A random parameter scan has been performed, 
which checks for each point, if all experimental constraints, 
including the relic density, are fulfilled. One observes that 
the selected region of Fig. [2] requires values of tan/3 above 
50, if the relic density constraint is to be fulfilled as well. 
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Fig. 4. flh 2 as a function of tan (3 is plotted for mo = 1500 
GeV, m 1 / 2 = 200 GeV, Ao = 0. The dependence on the 
uncertainty of the SM parameters (a s , m& and rat) covers 
values in a range of more than one order of magnitude. The 
horizontal shaded band is value of the Standard Model of 
cosmology. 



accordingly. For large values of tan/3, typically above 
50, this is the case, as shown in Fig. [3] In Fig. Q] the 
relic density is shown as a function of tan (3 for a small 
value of mi/2 which is compatible with the ERGET 
excess. The calculated value of the relic density has a 
strong dependence on the SM parameters like a s , m t 
and nib and the corresponding uncertainty covers val- 
ues in a region which spans over more than one order 
of magnitude. For large values of tan/3 the relic den- 
sity is very sensitive to tan /3 and the SM parameters, 
which can be understood in the following way: 



In mSUGRA models the parameters are set at the 
GUT scale. The quantities like cross sections or masses 
have to be calculated at much lower energies, typically 
the electroweak scale. Therefore renormalization group 
equations are used to run the couplings and soft break- 
ing terms from the unification scale to lower energies. 
This running also depends on the values of the SM pa- 
rameters in particular on the Yukawa couplings of the 
third generation and the strong coupling constant: 
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with t — \og(M 2 VT /Q 2 ) and the initial conditions mn 2 - 
mn 1 = rriQ . The Yukawa couplings h are determined by 
the running quark masses. For the running masses it is 
important to include higher order radiative corrections 
[8] and supersymmetric contributions [9], since they 
are not negligible. In particular corrections to mj, be- 
come more and more important with increasing tan/3. 
Since they are negative, the coupling is getting large 
and the RGE equation for mn 1 is dominated by the 
corresponding /i 2 -term. So a difference in the measured 
quark mass is assigned to the mass of the pseudoscalar 
Higgs mass by the definition of tua (at tree level): 
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Electroweak symmetry breaking is triggered if the fol- 
lowing condition is fullfillcd: 
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The parameter mi depends on the RGE of m# 2 . For 
large values of mo the running of to# 2 is very steep and 
sensitive to to* [5]. The reason for this is the large ht 
and the heavy sfermion masses, since Too is large, i.e. 
the term in the RGE of m# 2 proportional to h 2 gets 
dominant. If tan/3 is small than hb is also small and 
the running of mn 1 is not steep, since the term which 
includes the sfermion masses is small. In this case mjy 2 
turns out to have a value determined by the condition 
above with a relatively small theoretical uncertainty, 
since Mz, tan/3 and mi have small errors. 

In case of large tan (3 the running of vriHi starts 
to be steep too, due to the increasing bottom Yukawa 
coupling. The steep running of to# 2 and mn 1 has as 
an effect that from the minimization condition the pa- 
rameters to 2 and toi cannot be determined with high 
accuracy since two of the quantities in the relation 
have large uncertainties. This causes a high sensitiv- 
ity to the radiative corrections of the Higgs potential. 
Only a small change of these corrections or the heavy 
quark masses will result in very different values of toj^ 1 
and to# 2 and correspondingly via has large theoretical 
uncertainty. 

For such large values of tan (3 the annihilation via 
pseudoscalar Higgs exchange, being proportional to 
(tan/?) 2 , becomes dominant, if the neutralino mass is 
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Table 1. Typical mSUGRA parameters from the EGRET 
analysis, sparticle masses and observables. 
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not close to half the Z° mass. In the latter case the 
annihilation is almost a pure s-wave, i.e. at low ener- 
gies the cross section is suppressed. Therefore the fact 
that the annihilation is still observable in the present 
universe, as indicated by the EGRET excess, implies 
a lower limit of about 50 GeV on the neutralino mass, 
unless one allows for an extreme boosting of the DMA 
by the fact that the DM is clustered. One would need 
boost factors of several thousands [5] • 

A typical mass spectrum for the allowed region 
of Fig. and tan/3 = 52.2 is given in Table [Q It 
should be noted that the annihilation is still in the so- 
called bulk region, i.e the regions not dominated by co- 
annihilation or resonances, since the neutralino mass is 
far away from the pseudoscalar Higgs mass resonance 
for the mSUGRA spectrum and not close to any of 
the other sparticles, like stau or chargino, as shown 
in Table [T] The total annihilation rate, which is the 
sum of the self-annihilation and co-annihilation rate, 
is fixed by the observed relic density. Therefore a large 
co-annihilation rate automatically implies a negligible 
self-annihilation rate. Since in the present universe the 
NLSPs have decayed, only the self-annihilation is op- 
erative now and would be practically zero in case of 
strong co-annihilation. So it is good for indirect DM 
detection that the combination of EGRET data with 
Higgs mass limits results in a spectrum, for which the 
co-annihilation is negligible. 



3 Conclusion 

In our previous paper [2] the observed excess of dif- 
fuse Galactic gamma rays was shown to exhibit all the 
features of Dark Matter Annihilation, especially the 
spatial distribution of the excess was shown to trace 
the DM distribution, as proven by the fact that one 



could reconstruct the peculiar shape of the rotation 
curve of our Galaxy from the gamma ray excess. In this 
paper the DM interpretation of the EGRET excess is 
compared with Supersymmetry and it is shown that 
the minimal supersymmetric model with the popular 
supergravity inspired symmetry breaking, gauge uni- 
fication and radiative electroweak symmetry breaking 
is in perfect agreement with the EGRET excess. The 
mass spectrum of the gauginos is governed by the neu- 
tralino mass corresponding to iriiw roughly between 
125 and 175 GeV, while the scalar masses are con- 
strained by the Higgs mass and/or Br(B — > X s ~f) to 
have mo above 1200 GeV and to be below roughly 2 
TeV in order to allow for EWSB. The upper limit has 
a large uncertainty from the Yukawa coupling of the 
third generation and the gauge couplings, since the 
steep running of the mass terms of the Higgs potential 
are governed by them. Such a mass spectrum is observ- 
able at the LHC. If confirmed, especially a neutralino 
mass around 60 GeV, then this would prove that DM 
can indeed be considered the supersymmetric partner 
of the Cosmic Microwave Background, since the neu- 
tralino is almost a pure bino in this case. 
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